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Abstract 

We study the quasi-long-range ordered phase of a 2D XY model with quenched 
site-dilution using the spin-wave approximation and expansion in the parameter 
which characterizes the deviation from completely homogeneous dilution. The re- 
sults, obtained by keeping the terms up to the third order in the expansion, show 
good accordance with Monte Carlo data in a wide range of dilution concentrations 
far enough from the percolation threshold. We discuss different types of expansion. 
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The XY model in two dimensions is the simplest example of a system exhibit- 
ing "quasi-long-range order" (QLRO), which appears at low temperatures in a 
number of physical models of great importance, e.g. magnetic films with planar 
anisotropy, but also thin-film superfluids or superconductors, two-dimensional 
solids, 2d-classical Coulomb gas or fluctuating surfaces and the roughness 
transition [1,2]. Although no exact solution exists for this model, most of its 
properties are known from different approaches. 

Destruction of long range ordering is due to the presence of stable topologi- 
cal defects (vortices) [3,4,5], a situation which strongly contrasts with usual 
ordering in systems undergoing a ferromagnetic phase transition. First, the 
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magnetization of the XY model on an infinite 2D lattice remains zero at any 
non-zero temperature [6], thus it is impossible to describe in the thermody- 
namic limit the quasi-long-range ordered phase by this usual order parameter, 
however the spin-spin pair correlation function gives a distinct indication of 
QLRO. Its asymptotic behaviour changes from exponential at high tempera- 
tures to power low decay at low temperatures. This transition is referred to 
as the Berezinskii-Kosterlitz-Thouless (BKT) transition and the point where 
this change of behaviour occurs is the BKT temperature. 

A quantity of interest which characterizes the QLRO phase, is then the tem- 
perature dependent exponent of the correlation function: 

rj(T) = - lim ln(Sr ^ +R) . (1) 
/v ; |R|-oo In R V ; 



The spin-wave analysis of the model gives a reliable value of i] for small enough 
temperatures [8]. The reliability of the harmonic approximation for this model 
is grounded by the RG analysis [5]. 

The case of the (classical) XY model on a regular lattice (without defects) 
ff = -^EE^- O (5?5* + S*5* ) , (2) 



has been studied intensively (when J(r — r') is limited to nearest neighbours) 
and its properties are well known (see e.g. Ref. [2]). The addition of defects 
(site-dilution, bond-dilution) has been considered as a trivial modification, 
since the Harris criterion [9] states in this case that the universality class of the 
diluted model remains the same as that of the pure one. It means that the crit- 
ical exponents of both pure and disordered models are unchanged, when eval- 
uated at their corresponding BKT points, e.g. r/ dil (T^ T ) = V 1 (Xbkt), but 
the functions rj pure (T) and 77 dll (T) characterizing the low temperature phase of 
pure and disordered systems are different and the exact behaviour of ?7 dll (T) 
is a question which deserves attention. For example, it is not obvious how 
the impurities can interact with the vortices and influence the QLRO. This 
question is addressed e.g. in Refs. [10,11]. 

In a previous paper [12] the influence of uncorrelated (normally distributed) 
site-dilution was considered and a perturbation expansion for the case of weak 
dilution was proposed. The two-spin coupling term J(r — r') in Eq. (2) was 
replaced by J(r — r')c r c r / with 

1, if the site r has a spin: 

(3) 

0, if the site r is empty. 



2 



The physical quantities which characterize the system with quenched disor- 
der after the thermodynamical averaging must be averaged over all possible 
configurations of dilution. This configurational averaging is denoted as (...): 

ro = n e [^i- Cr ,o+(i-c)5 Cr , ](...), (4) 

r c r =0,l 

where c is the concentration of occupied sites. 

Starting with the Hamiltonian in the harmonic approximation, 

H = 7 E E J ( r " r ') (*r " Or'? c r c r , , (5) 

4 r r ' 

and realizing the Fourier transformation of the variables: 

^• = 4=E ^ = ^=E e~ lkr °r, (6) 

ViV k ViV r 

j(r) = iE^Kq), Kq) = E^W. ( 7 ) 

q r 

(N is the number of sites in the lattice, and k runs over the 1st Brillouin 
zone), one has 

H = JJ2 TkM-k + J E E(7k+k' - 7k - Tk'Mk + k')0 k 0k' (8) 

k k k' 

+ J E E E(2 - 7q) [P("k - k' - q)p(q) - p(-k - q)p(-k' + q)] Mk' 
k k' q 

where p(q) = X)r e_iqr (l — c r), and 7 k = ^7 [^(0) — ^(k)] = 2 — cos k x a — 
cos k y a on the square lattice. The first term is the Hamiltonian of the pure 
system, so one can consider p as a parameter of perturbation of this Hamil- 
tonian. Note, that power of p corresponds to the number of sums over k. A 
classification of the perturbation theory series with respect to number of sums 
over k corresponds to the expansion in the ratio of the volume of effective 
interaction to the elementary cell volume [13]. Taking this ratio to be small 
means that it is valid for the short-range interacting systems, which holds for 
our problem. 

The linear approximation in p-expansion presented in Ref. [12] gives the result 
for the exponent of the pair correlation function: 

V dil = (1 + 2(1 - c)) r] purc + 0{p 2 ). (9) 
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Here we can report the result for this expansion up to the second order in p: 
r] dil = (l + 2.73(1 -c) + 1.27(1 - c) 2 ) ?? purc + 0(p 3 ), (10) 

The figures follow from numerical estimate of sums. The 1st and 2nd order 
perturbation expressions fit the Monte Carlo results only for very small con- 
centrations of dilution (see Fig. 1). Unfortunately the calculation of the third 
order seems to be too tedious and perhaps does not deserve so much effort, 
therefore it is desirable to investigate another road. 

In the present paper we propose to introduce the parameter of expansion in a 
different manner in order to extend the region of reliability of the expansion to 
stronger dilutions (but of course still far enough from the percolation threshold 
where the whole approach fails) and to improve convergence. We will keep the 
notation p, but from now on one should understand it as the deviation from 
homogeneously diluted system: 

P(q) = ^E^ qr (Cr- C ). (11) 



We did not make any assumption about weakness of disorder, we may thus 
expect that the results of this expansion will be less sensitive to the value of 
dilution c. Rewriting the Hamiltonian with this new parameter one gets 

# = c 2 J ]T 7kM-k - cJ ]T ]T(7k+k' - 7k - 7k ,)p(k + k')0 k 0k' 

k k k' 

+ iEEE( 2 - 7q) [p(-k - k' - q)p(q) - p(-k - q)p(-k' + q)] k k , 
k k' q 

= c 2 H purc + H p + H p 2 , (12) 

where the first term is the Hamiltonian of the pure system now with a renor- 
malized coupling. 

The spin-spin pair correlation function, 



cos(# r+R - 6 r )e-P (Hp+H p 2) 



G 2 {R) = c r c r+R (cos(6> r+R - P )) = c r c r+R / -p(h p +h 2 )\ * ( 13 ) 



can be expanded in p and then configurationally averaged. In (13), (...)„, 
stands for the averaging with the pure system Hamiltonian with renormal- 
ized coupling c 2 J. Using (4) one has the equalities: p(q) = 0, p(q)p(q') = 
c{\ - c)i5 qWi0 and p(q)p(q')p(q") = c(l - 3c + 2c 2 ) ^ 5 q+ q'+q",o, which, 
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inserted into the expansion lead after tedious calculations to the third order 
expression 



c r c r+R (cos(6» r+R - r )) = (? (cos(6> r+R - r )) 

X 



, _ 1 - c 1 ^ (7k+k' - 7k ~ 7k') 2 . 2 ^ 
c 3 PJN*^, 7 2 7k , Sm 2 



l-3c + 2c 2 1 / 



_2_ „ sm — 

/3J\ N ^ k ^ 

(7 k+k ' -7k-7 k 0(7k' -k" -Tk' ~7k" )(7 k +k" -7k~7 k ") 2 kR 

~ ~ ~" ~ ' " " bill . \ 



N 3 ^k,k',k" 7k 7k' 7k" 

In the limit iV — > oo, R — > oo we have found for the sums in (14): 



(14) 



J_ — 2_ pj, cons t _)_ J-\ n E 

AT ^ 7k 2?r a ' 



k 



-^VV (7k+k '7 k " v)2 sin 2 ^ » const' + ^ln*, (15) 

yy2 Z-^ Z-^ 7kTk' 2 2?r a 

k k' 

J_ (7k+k'-7k-7k')(7k'-k"-7k'-7k")(7k+k"-7k-7 k ") - 2 kR 

N 3 , ?Tt „ ^7k'7k" 2 



k,k',k" 



« const" + S^Zln^. 

The figures 0.73 and 0.27 come from numerical summation. It appears that to 
zeroth-order, the change of exponent comes from a renormalization of the cou- 
pling strength, the first-order term is identically vanishing. For small enough 
temperatures it is now possible to write the pair correlation function in the 
power law form: 

G 2 (R) « c 2 (i?/a)-" dU . (16) 



Reminding the correlation function exponent of the pure system in the SW 
approximation, i] puTC = (27T/3J) -1 , we can write 

^ dil = ?T rc + °- 73 ^ - 0-27 1 ~^ 4 +2c2 ) + 0(p 4 ). (17) 

The first term in the brackets, 1/c 2 , corresponds to the zeroth order in the 
expansion, the first-order term is identically vanishing as was already noted 
before, the second and third terms in the brackets correspond to the second- 
and third-order terms in the p-expansion respectively. 
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Fig. 1. Comparison between the 1st (which in fact coincides with the 0th), 2nd 
and 3rd-order expansions and MC simulations at very low temperatures (the values 
of {(5J)~ l are indicated in the legend). Expressions (9) and (10) of the previous 
expansion are shown (thin lines) for comparison. Insert shows the vicinity of the 
pure system. 

In order to check these expressions, simulations of 2D XY-spins are performed 
using Wolff's cluster Monte Carlo algorithm [15]. The low-temperature phase 
being critical, local updates of single spins would suffer from the critical slow- 
ing down. Implemented in the case of the XY model, the Wolff algorithm first 
introduces bonds through the Ising variables denned by the sign of the projec- 
tion of the spin variables along some random direction. Then clusters of sites 
are built by a bond percolation process (here the random graph model of the 
Fortuin-Kasteleyn representation). The percolation threshold for these bonds 
coincides with the Kosterlitz-Thouless point [16], which guarantees the effi- 
ciency of the Wolff cluster updating scheme [15] at Tkt- In the low-temperature 
phase we are interested in, this algorithm could be less efficient, but neverthe- 
less preferable to a local updating, since the correlation length is diverging. 
Using this procedure, we discard typically 10 5 sweeps for thermalization, and 
the measurements are performed on typically 10 5 production sweeps. Averages 
over disorder are performed using typically 10 3 samples. There is no need of a 
better statistics. The boundary conditions are chosen periodic and the critical 
exponent rj(T) of the correlation function is measured indirectly through the 
finite-size scaling behaviour of the magnetization 

M T {L)~L-^ T \ x a (T) = ^(T), (18) 
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where the last scaling relations holds in two dimensions. 

In Figure 1, we compare the Oth to 3rd order expansions (remember that 
the lrd order term vanishes) with the MC data. The agreement is quite good 
using the expansion parameter (11) which provides a clear improvement of the 
previous expansion given by expressions (9) and (10). Of course the question 
of the next order is not settled, but now we are on the way to counting higher 
orders or even summing the whole series. Another direction of future work 
would be to implement the same type of perturbation expansion within the 
Villain model [17] and to explore the deconfming transition of the diluted 
model. 
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